Modulated superconductivity near Pomeranchuk instabilities in the 

spin channel 
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We study the competition between a Pomeranchuk instability in the spin channel with angular 
momentum I — 1 and an attractive interaction, favoring Cooper pair formation. We find, at 
mean-field approximation, that superconductivity strongly suppress the Pomeranchuk instability. 
Moreover, we have found a metastable modulated superconducting phase with similar characteristics 
of the FFLO state. 
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A Fermi liquid is, except in one dimension, a 
very stable state of matter—. At least two types of 
instabilities, related with attractive interactions, are 
known: Pomeranchuk 2 and superconducting instabilities. 
Pomeranchuk instabilities occur in the presence of two- 
body interactions containing a strong attractive compo- 
nent in the forward scattering channel with a definite 
angular momentum. In the context of Landau theory, 
the instability sets in when one or more dimensionless 
Landau parameters with angular momentum £ in 

the charge (s) or spin (a) channel, acquire large negative 
values. 

Pomeranchuk instabilities in a charge sector sponta- 
neously break rotational symmetry. In particular, an in- 
stability in the F 2 S channel produces an ellipsoidal defor- 
mation of the Fermi surface^. From a dynamical point 
of view, the resulting anisotropic ground state is a non- 
Fermi liquid, due to overdamped Goldstone modes that 
wipe out quasi-particle excitations^— . 

In the spin channel, the ferromagnetic Stoner 
instability 8 occurs when Fq acquires large negative val- 
ues, producing a divergence in the magnetic susceptibil- 
ity. This phase transition preserves rotational symmetry, 
however it breaks time-reversal symmetry. Higher order 
angular momentum interactions produce anisotropic as 
well as isotropicphases. Several examples were studied 
in detail in refs. 0and[T3- The 1=1 channel have special 
interest. When F" < 0, an ordered isotropic and time- 
reversal invariant phase is possible. This phase, called 
/3-phase in ref. ; &, dynamically generates a spin-orbit 
coupling. This is a very interesting possibility, since it 
allows the generation of spin-orbit couplings from many- 
body correlations, differently from the usual one-particle 
relativistic effect. 

On the other hand, superconductivity is developed 
in the presence of a small attractive interaction in the 
particle-particle (BCS) channel. The superconducting 
(SC) state is generally characterized by a complex order 
parameter which breaks gauge symmetry, A at0 -> (r, f") = 
{^(r)^,^)), where the operator ip^{f) creates an elec- 
tron with spin a at the position r. The usual classi- 



fication of Ag-.o-' (r, f) as s-wave, d-wave, p-wave, etc, 
resides in the irreducible representations of the lattice 
point group. Also, the absence of spin-orbit interactions 
allows the additional differentiation between singlet and 
triplet order parameters. However, the SC state could 
break lattice translation and/or rotational symmetry. In 
that case, this classification is no more possible. One 
particular example is an oscillating order parameter like 
Ao- !(T ' (r, 0) = A CTjCr / cos(<f • r) , where q is an ordering wave- 
vector. In a modulated superconducting state, proposed 
by Fulde, Ferrell, Larkin and Ovchinniko v 11 ! 12 (FFLO 
state), the spatial modulation of the Cooper pairing is 
due to a mismatch of Fermi surfaces, produced by an ex- 
ternal magnetic field (by Zeeman effect). FFLO states 
have also been proposed to occur in other scenarios like 
in imbalanced cold atoms with different species^, and in 
heavy fermions systems when different orbitals hybridize 
under external pressure^. Also, for generally non-local 
BCS potentials, a modulated SC order parameter can 
coexist with charge density waves^. Recently, a striped 
order parameter, called "Pair Density Wave, (PDW)", 
was proposed^— to explain anomalous transport prop- 
erties observed in cuprates superconductors. 

Although Pomeranchuk and BCS instabilities are gen- 
erated by attractive interactions, they are competing 
ones. A superconducting gap suppresses Fermi surface 
deformation. A detailed example of this effect was stud- 
ied in refs. [Til and [2(| where the competition between 
d-wave Pomeranchuk instabilities and d-wave supercon- 
ducting order parameter was considered. 

In the spin channel, a stronger competition is expected, 
since, in general, magnetic order is expelled from a su- 
perconductor state. However, it is possible to have in- 
stabilities with higher angular momentum that preserve 
time-reversal invariance. In this paper, we analyze an 
example of this class of systems. In particular, the (3- 
phase introduced in ref. [9|, opens the possibility of the 
formation of Cooper pairs with zero helicity and finite 
momentum, possibly producing a modulated supercon- 
ducting state. 

The main point of this paper is to report on the com- 
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FIG. 1. Phase diagram. Ao is the dimensionless isotropic 
BCS gap A = bA /v F k F . Ff = /fiV(O) is the dimensionless 
antisymmetric Landau parameter with angular momentum 
I = 1. The solid line, separating the /? and the SC phase, 
is a first order line. The dashed line indicate the onset of a 
metastable modulated SC phase. 



peting character of Cooper pairing and a Pomeranchuk 
instability in the I = 1 spin channel. Our model is 
based on a 2D Fermi liquid with interactions in the 
Fy channel and a general attracting interaction favoring 
Cooper pairing. The main results are summarized in the 
phase diagram of fig. (JT]), obtained in a self-consistent 
mean-field approximation described below. In the figure, 
Ao ~ e" 1 ^ 9 is a parametrization of the BCS coupling con- 
stant g. Physically, it represents the gap of the uniform s- 
wave superconductor in the absence of any Pomeranchuk 
instability. For Ff > —2, the system is in the usual 
s-wave SC state. When F® < —2 the Pomeranchuk in- 
stability sets in and the mean-field ground state at low 
Ao is the so called /3-phase consisting on a Fermi liquid 
with an induced spin-orbit interaction. For higher values 
of Ao, s-wave superconductivity strongly suppresses the 
/3-phase. They are separated by a coexisting line, signal- 
ing a first order phase transition. Moreover, for greater 
values of the BCS coupling, a metastable modulated su- 
perconducting phase appears with a finite wave- vector q, 
given roughly by \q\ ~ kp\/l — 2/\F{ a \. 

In the rest of the paper we present our model and 
sketch the mean-field calculation leading to the phase 
diagram of fig. (p}. A spin-current operator with orbital 
angular momentum I = 1 is given by 



V af *(x) = {x)a^V a ^(x) 



(1) 



where ip is a two-component spinor, a 11 are the Pauli 
matrices, and a = x,y. Ve lkx = i^e lkx . 

Following ref. we write a Hamiltonian with forward 



two-body interactions in this channel 

H = [ d 2 x ip^x) (e(-iV) - (i) ip(x) + 



+ -J d 2 xd 2 x'f1{x - x')V a »{x)V afi (x') . (2) 

The Fourier transform of /f(x) is given by /f(fc) 
" 4- n\f^\k 2 ), defining, in this way, an effective in- 
teraction range r = yftf/f]. For a particle-hole sym- 
metric system we consider the following expansion of the 
dispersion relation around a circular Fermi surface, 



e(fc) — fi = v F ■ [k — kp] ■ 



(v F k F y 



v F -[k- k F ] ) , (3) 



where the dimensionless parameter b measures the effec- 
tive curvature of the band near the Fermi surface and we 
have ignored terms proportional to (k — k F ) 5 . 

The /3-phase is defined by the mean-field Hamiltonian 



Hmf — 



d 2 k 

WT 2 



^(k)(e(k)-n-na-k)ij(k), (4) 



where n is determined self-consistently by 



1 



d 2 k 
(2tt)2 



(^(k)(a-k)^(k)) . (5) 



This mean- field theory is valid when k F y/n\fi \ ^> 1, i. 
e. , when the range of the interaction is much larger than 
the interparticle distance. 

In references d and [l(| eq. ([5]) was solved by con- 
sidering a Fermi liquid ground state. When n / the 
spectrum of Hmf splits into two opposite chiralities with 
dispersions = e(k) — (fi + n) and = e(k) — {p, — n), 
as shown in figure ([2]). The Fermi momentum of each 
branch is given by kp = kf ±g/2 and the relation be- 
tween q and n is computed using eq. ([3]), 



2 



(v F k F y 



0(n 5 ) 



(6) 



(we have renormalized the chemical potential to keep the 
density constant). 

The ordered phase is characterized by a spontaneously 
generated spin-orbit coupling with the global rotation in- 
variance unbroken. The spin and orbital angular momen- 
tum are not conserved independently, however the total 
angular momentum J — L + S is conserved. In the he- 
licity base C = (Ct; Cj,)i where the operator k ■ a is diago- 
nalized with eigenvalues ±1, the mean- field Hamiltonian 
takes the simpler form 



Hmf — 



d 2 k 



(2tt)2 

with the self-consistent equation 



CHk){e{k)-^-na z )C{k), (7) 



(8) 
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FIG. 2. Fermi Surface splitting in the helicity base and pairing 
in the /3-phase. 



A BCS type interaction in the ip, ip* base induces differ- 
ent types of interactions in the /3-phase, where the usual 
Fermi liquid properties are lost. In the helicity base (, C , 
intra-band p-wave pairing, as well as inter-band finite 
momentum pairing, are induced. These interactions are 
competing with each other and with interactions in the 
forward scattering (Pomeranchuk) channel. In this work, 
we study the competition between an inter-band SC or- 
der parameter and the /3-instability, leaving a complete 
study of the phase diagram, including other Cooper pair 
symmetries and fluctuations, for a future presentation. 

We consider a simple model where we introduce the 
following attractive interaction favoring Cooper pair for- 
mation with linear momentum q, 

d 2 k 



H 



sc 



Cl.-j 



t,fc+qr/2H,-fc+g/2 



h.c. 



(9) 



Note that q (eq. ©) is given by the Pomeranchuk in- 
stability and should be determined sclf-consistently. The 
"gap equation" that complements eq. © is 



A, = -g 



d 2 k 



(<4,-fc+9/2(t,fc+«/2) 



(10) 



(2tt) 2 

To find a self-consistent solution, it is necessary to solve 
eqs. ©, (|10l) and © simultaneously, where the ex- 
pectation values are taken with the Hamiltonian H = 
Hmf + Hsc- 

Diagonalizing H, considering /" < 0, we find at zero 
temperature 

m, / ^ 

n 



4 / j^M-^)-^-)} , (11) 



A 9 = -g 



d 2 k A n 



{e(-w+) - 0(-w_)} 



(12) 



(2tt)2 2£ 

where O is the usual Heaviside function and single- 
particle excitations read 



w± 



' t _ 4- 

£ k+q/2 6 k-q/2 



\ 



' t , I 

6 k+q/2 > E k- 



q/2 



\W 2 



(14) 



It is simple to realize that the s-wave BCS order param- 
eter n = 0, A g = A is always a solution of eqs. (ITU and 
(fT2"l) . On the other hand, in the absence of superconduc- 
tivity, Ao = 0, a solution with finite n ~ q has lower en- 
ergy for F£ < —2, leading to the /3-phase. In this regime 
and for small values of the BCS coupling constant, the j3- 
phase has always smaller energy than the uniform super- 
conducting one. For stronger superconducting coupling, 
the /3-phase is suppressed and the mean-field ground 
state is a uniform superconducting phase, as shown in 
figure ([TJ. The first order line that separates these two 
phases was evaluated by equating the mean-field Hamil- 
tonians, (H)a = (H)p, computed in the two possible so- 
lutions: uniform superconductor {n = 0, A q = Aq} and 
the /3-phase {n ~ (v F k F /b 1 / 2 )^/! - 2/\Ff\, A q = 0}, re- 
spectively. Note that the band curvature at the Fermi 
surface (6) is essential to stabilize the /3-phase. 

Now, we explore the interesting possibility of a mod- 
ulated solution n^g^OjAg^O. The fc-integrals in 
equations (fTTj) and (|12[) are strongly constrained by the 
Heaviside functions. The main contribution to eq. (|11[) 
comes from the region w + < 0. Written in polar coordi- 
nates (k,9), with cos 9 = {k ■ q)/kq, and first integrating 
over k, we get 

n = \F?\n 

.a- 2ir 




(v F k F ) 



0) (1 + 



cost 



(15) 



where we have defined the dimensionless quantity x = 
A q /n. 9 (x, n) are the integration limits which keep the 
argument of the square root positive. Near the criti- 
cal point, we can expand the integrand in powers of 
n/v F k F « 1, obtaining 



- Fi - 
n = ——n 



C (x) 



[vf kf 



Tn 2 d(x) 



0((- 



-) 5 ), 



y v F k F 

(16) 

with C (:e) - 1 - (x/2) 3 / 2 , C\(x) = 1 + x ln(l + x/S) and 
Fi = /i N(0), where N(0) is the density of states at the 
Fermi surface. 

On the other hand, the integral in eq. (TT21 has the 
usual ultraviolet divergence of the BCS gap equation. A 
convenient way to deal with this integral is to sum and 
subtract 0(u;_) in the second term and then to subtract 
the identity 



1 



9 



d 2 k 1 
(2tt) 2 2£° 







(17) 



(13) from the first term of eq. (TT21 . £° = ^/e 2 + \A \ 2 and 
Aq is the uniform superconducting gap in the absence of 
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Fermi surface splitting. With this, eq. (|12[) is rewritten 

as 



AT(O) 



In 



A 



d 2 k 1 
(2^2£ 



{0(-w+) + 6(u;_)} . (18) 



In this way, the above integral is controlled. The cou- 
pling constant g and the ultraviolet cut-off are both con- 
tained in the definition of Ao (eq. (fl7j)). As before, 
near the Fermi surface, where interactions are important, 
0(cj_) = 0. Performing the remaining integral over k, 
the gap equation is written as 



In 



A n 



A, 



r(x) 



(19) 



where x = A q /n, and we have defined the function 



T(x) 



S(x) 



2tt 



sinh 



cos ( 



(20) 

with Q(x) = cos^ 1 (a; - 1) for x < 2 and T(x > 2) = 0. 
Therefore, to look for modulated superconductivity we 
need to solve eqs. (HI))) and (|L9")) self-consistently. Firstly, 
we note that for x > 2 the only solution, {n — 0, A q — 
Ao}, corresponds to the uniform superconducting phase. 
This imposes a lower limit for modulated superconduc- 
tivity since, for obtaining a non-trivial solution, we must 
have A 9 < 2h < Aq. However, for any value of the pa- 
rameter Ao > 2n, the uniform superconducting mean- 
field energy is always lower than the modulated one, 
{H)a < (H)a ■ On the other hand, there is a region of 
the phase diagram, shown over the dashed line in figure 
([1]), in which (H)/\ q < (H)p. Therefore, modulated su- 
perconductivity appears as a metastable phase for finite 
values of Ao and \Fi\ > 2. The precise location of the 
onset of metastability depends on the relation between 



the band-width, in which interactions are relevant (the 
energy cut-off), and the band-curvature at the Fermi sur- 
face. 

Summarizing, we have studied, in the context of Fermi 
liquids, the competition between a Pomeranchuk insta- 
bility in the spin channel with angular momentum I = 1, 
triggered by the Landau Parameter and an attractive 
interaction favoring Cooper pair formation. We built up 
a zero temperature mean-field phase diagram in terms of 
the two parameters of the model: F® and g, parametrized 
by Ao. We found that, as observed in other channels 
of Pomeranchuk instabilitie s) 19 ' 20 , the superconducting 
gap reduces the tendency to Fermi surface deformation. 
Moreover, we have found a metastable modulated super- 
conductor, enhanced for larger values of \F"\. 

The phase diagram of figure ([!} could be modified by 
several reasons. On the one hand, we have ignored intra- 
band homogeneous pairing that surely competes with the 
inter-band pairing considered in this paper. On the other 
hand, long wave-length fluctuations will certainly modify 
mean-field results. There are two types of important fluc- 
tuations to be considered. In the /3-phase, there are three 
branches of Goldstone modes transforming with SO(3) 
group, associated with the breaking of relative spin-orbit 
symmetry^. The longitudinal mode has a linear disper- 
sion relation while the transverse ones are responsible for 
a Lifshitz instability, signaling a tendency to a modulated 
ground-state with parity symmetry breaking. This effect, 
in the presence of pairing interactions, should probably 
enhance the modulated SC order. 
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